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We study several cosmological models with Bianchi VIo symmetries under the self-similar ap- 
proach. In order to study how the "constants" G and A may vary, we propose three scenarios where 
such constants are considered as time functions. The first model is a perfect fluid. We find that 
the behavior of G and A are related. If G behaves as a growing time function then A is a positive 
decreasing time function but if G is decreasing then A is negative. For this model we have found 
a new solution. The second model is a scalar field, where in a phenomenological way, we consider 
a modification of the Klein-Gordon equation in order to take into account the variation of G. Our 
j third scenario is a scalar-tensor model. We find three solutions for this models where G is growing, 

„ , constant or decreasing and A is a positive decreasing function or vanishes. We put special emphasis 

on calculating the curvature invariants in order to see if the solutions isotropize. 
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I. INTRODUCTION 



Current observations of the large scale cosmic microwave background (CMB) suggest to us that our physical universe 
is expanding in an accelerated way, isotropic and homogeneous models with a positive cosmological constant. The 
analysis of CMB fluctuations could confirm this picture. But other analyses reveal some inconsistencies. Analysis of 
WMAP data sets shows us that the universe might have a preferred direction. For this reason, it may be interesting 
to study Bianchi models since these models may describe such anisotropics. 

The observed location of the first acoustic peak of the temperature fluctuations on the CMB corroborated by the 
data obtained in different experiments [l[ , indicates that the universe is dominated by an unidentified "dark energy" 
' and suggests that this unidentified dark energy has a negative pressure . This last characteristic of the dark energy 
points to the vacuum energy or cosmological constant as a possible candidate for dark energy. 

In modern cosmological theories, the cosmological constant remains a focal point of interest (see 0-@ for reviews 
of the problem). A wide range of observations now compellingly suggest that the universe possesses a non-zero 
cosmological constant. Some of the recent discussions on the cosmological constant "problem" and on cosmology 
with a time- varying cosmological constant point out that in the absence of any interaction with matter or radiation, 
■ the cosmological constant remains a "constant" . However, in the presence of interactions with matter or radiation, a 
solution of Einstein equations and the assumed equation of covariant conservation of stress-energy with a time- varying 
A can be found. This entails that energy has to be conserved by a decrease in the energy density of the vacuum 
component followed by a corresponding increase in the energy density of matter or radiation. Recent observations 
strongly favour a significant and a positive value of A with magnitude A(Gfi/c 3 ) ps 10~ 123 . These observations suggest 
an accelerating expansion of the universe, q < 0, @- 

Following Maia, et al 7] who have pointed out that although the cosmological A-term has a very small value today, 
it may contribute to the total energy density of the universe. For this reason, since its present value, Ao, may be a 
remnant of a primordial inflationary stage, it seems natural to study cosmological scenarios which include a decaying 
vacuum energy density in such a way that it must be high enough at very early times and sufficiently small at present 
times in order to be compatible with the current observations. One of the first attempts at considering a decreasing 
cosmological term was formulated by Chen et al Q . By studying the Wheeler-DeWitt equation, they argue through 
dimensional considerations that the cosmological A-term must follow a relationship such as A ~ t~ 2 , in order to fit 
with current observations. Other mechanism to describe such variation have been formulated within the framework 
of the so-called "quintessence models" . Recently this class of models have received a great deal of attention Q and 
[Toj . Taking into account different observational data it is possible to rule out and to obtain "correct" potential which 
could play the role of an effective cosmological constant. This strengthens the idea of considering alternative theories 
where the scalar field is non-minimally coupled to gravity, like scalar-tensor theories (STT) [llj . This class of theories 
furthermore allows the variation of other constants such as the Newton gravitational one. There are several STT 
derived from the original one, the Brans-Dicke (BD) model (see for example [Hj]-[l7j])- They have been formulated 
as possible solutions to the discrepancies with observations and try to explain the behaviour of the universe at late 
times (see [ijj]- [HJ). Of particular interest are the so called chameleon scalar-tensor theories (22j . 

The study of self-similar (SS) models is quite important since a large class of orthogonal spatially homogeneous 
models are asymptotically self-similar at the initial singularity and are approximated by exact perfect fluid or vacuum 
self-similar power-law models. Exact self-similar power-law models can also approximate general Bianchi models at 
intermediate stages of their evolution. This last point is of particular importance in relating Bianchi models to the 
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real Universe. At the same time, self-similar solutions can describe the behaviour of Bianchi models at late times i.e. 
as t — y oo (see [23jj ) . 

The aim of this work is to study self-similar solutions of a Bianchi VIo cosmological model in different contexts 
and where the "constants" G and A may vary. We are mainly interested in finding exact solutions for the proposed 
models as well as to compare the behavior of G and A in the different contexts. In section 2 we start showing all 
the geometrical ingredients that we are going to use throughout the paper. We put special emphasis on the study 
of the curvature invariants in order to study whether the obtained solutions isotropize. Once we have calculated the 
homothetic vector field (HVF) in section 3 we study the "classical" solution for a vacuum and perfect fluid models 
comparing these solutions with those obtained ones in a previous work (in that work we used another Bianchi VIo 
metric [24j ) as well as a perfect fluid model with time-varying constants. In section 4, we start by studying the kind 
of potential and scalar fields compatible with self-similar solution. The stated theorems are very general and are 
valid for all the Bianchi models as well as for the FRW models. All the proofs have been performed by studying the 
Klein-Gordon equation through the Lie group method. Once we know the scalar fields compatible with the self-similar 
solution we continue studying a simple scalar model as well as a non-interacting scalar model with matter. In order 
to incorporate a gravitational "varying-constant" G(t) within this framework we purpose, in a phenomenological way, 
a modified Klein-Gordon equation. As above we need to study the class of potential compatible with a self-similar 
solution and a varying G. Two kinds of models are studied. In section 5 we study a generalized scalar-tensor model 
that determine an accelerated expansion at the present epoch, with arbitrary cu((f>) — const, and A((/>), where this last 
function plays the role of an effective cosmological constant. We would like to emphasize that in order to study the 
resulting field equations (FE) we have not needed to make any assumption, otherwise, we have deduced, the form of 
A(</>) by studying the conservation equation through the Lie group method. In section 6, we summarize our results. 
Finally, in the appendix A, we study through the matter collineation method the kinds of potentials compatible with 
a self-similar solution in the framework of G constant. In appendix B we study, using the same method, the G— var 
framework in such a way that we regain through this method the results obtained in section 4. 

II. THE GEOMETRIC INGREDIENTS 

We start by considering the following Killing vector fields (KFV) (see 26]) 

£i = d x + mzd y + myd z , £ 2 = d y , £ 3 = d z , (1) 

then 

Ei,6] = -Tn£ 3 , [&,&] =0, =m&. (2) 

Note that in this approach is essential to consider the m— parameter, otherwise it is impossible to obtain self-similar 
(SS) solutions. 

In this way it is obtained the following vector fields {Xj}, such that, Xj] = 0, = —C^Xk ■ 

X\ = coshmxdy + sinh mxd z , Xi = sinh mxd y + cosh mxd y , A3 = d x , (3) 

and the dual 1-forms: 

ui 1 = dx, (j 2 = coshmxdy — sinh mxdz, a; 3 = — sinh mxdy + cosh mxdz, (4) 
The metric is defined by 

ds 2 = -c 2 dt 2 + a 2 {t) (u l f + b 2 {t) (u 2 f + d 2 {t) (w 3 ) 2 , (5) 

finding that following metric when using Eq. (U) 

ds 2 = — dt 2 + a 2 dx 2 + (b 2 cosh 2 mx -I- d 2 sinh 2 mx) dy 2 

— 2 (b 2 + d 2 ) cosh mx sinh mxdydz + (b 2 sinh 2 mx -I- d 2 cosh 2 mx) dz 2 , (6) 

where we have set c = 1 . 

We may define the four velocity as follows: u % — (1, 0, 0, 0) , in such a way that it is verified, g(u z , u l ) = — 1. From 
the definition of the 4— velocity we find that: 

TT I (a 1 V d'\ 1 v^ rr d ( 3 \ , 
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and 

Isotropization means, in essence, that at large physical times t, when the volume factor, v — abd, tends to infinity, 
the three scale factors (a,b,d) grow at the same rate [27j • We will therefore say, by definition, that a model is 
isotropizing if, for each scale factors, a/ f — » const > 0, b/ f — > const > and d/ f — > const > 0, as v — > oo, where, 
/ = v 1 / 3 , is the mean scale factor. Then, by rescaling some of the coordinates, we can make a/f — > 1, bj / — > 1, 
d/f — > 1 and the metric will become manifestly isotropic at large t. Two such criteria are A — > and a — > 0, where 
the mean anisotropy parameter A is defined for the metric as (see, e.g., [28]) 



A ~3H 2 ~ 3^H* L (9) 

The mean anisotropy parameter gives a dimensionless measure of the anisotropy in the Hubble flow by comparing 
the shear scalar a to the overall rate of expansion as described by H . The anisotropy in the temperature of the CMBR 
enables one to estimate the value of a 2 at the present epoch. 

We also study the curvature behaviour of the solutions (see for example [29U3ll ] and HI])- The studied curvature 

3^0 i 

I4 = EijE l i , [33jj and the magnetic scalar I5 — HijH 1 ^ , of the Weyl tensor. The Weyl parameter W [33j, which is a 
dimensionless measure of the Weyl curvature tensor, 



quantities are the following ones: Ricci scalar, I = i?*, Krcstchmann scalar, Ii — R^mR 1 ^ 1 , the full contraction 
of the Ricci tensor, I 2 — RijR^ . The Weyl scalar, 7 3 = C ahcd C a bcd = h — 2/2 + as well as the electric scalar 



TI/2 1 T 

w2 =m = wi^ Eij + H ^ = 24W- (10) 

can be regarded as describing the intrinsic anisotropy in the gravitational field [34[. Cosmological observations can, 
in principle, give an upper bound on VV, although obtaining a strong bound is beyond the reach of present-day 
observations. 

Finally, we shall calculate the gravitational entropy. From a thermodynamic point of view there is every indication 
that the entropy of the universe is increasing. Increasing gravitational entropy would naturally be reflected by 
increasing local anisotropy, and the Weyl tensor reflects this. One suggestion in this connection was Penrose's 
formulation of what is called the Weyl curvature conjecture (WCC) [35| . The hypothesis is motivated by the need 
for a low entropy constraint on the initial state of the universe when the matter content was in thermal equilibrium. 
Penrose has argued that the low entropy constraint follows from the existence of the second law of thermodynamics, 
and that the low entropy in the gravitational field is tied to constraints on the Weyl curvature. Wainwright and 
Anderson [36j express this conjecture in terms of the ratio of the Weyl and the Ricci curvature invariants, 

p2 = T 2 (ID 

The physical content of the conjecture is that the initial state of the universe is homogeneous and isotropic. As 
pointed out by Rothman and Anninos [13, HH (see als o [39l0 the entities P 2 and I 3 are "local" entities in contrast 
to what we usually think of entropy. Gr0n and Hervik ([30|, [3l|) have introduced a non-local quantity which shows a 
more promising behaviour concerning the WCC. This quantity is also constructed in terms of the Weyl tensor, and it 
has therefore a direct connection with the Weyl curvature tensor but in a "non-local form". 

For SS spacetimes, Pelavas and Lake ([13]) have pointed out the idea that Eq. (fTTjl is not an acceptable candidate for 
gravitational entropy along the homothetic trajectories of any self-similar spacetime. Nor indeed is any "dimensionless" 
scalar. It is showed that the Lie derivative of any "dimensionless" scalar along a homothetic vector field (HVF) is zero, 
and concluded that such functions are not acceptable candidates for the gravitational entropy. Nevertheless [4l[ , since 
self-similar spacetimes represent asymptotic equilibrium states (since they describe the asymptotic properties of more 
general models), and the result P 2 — const., is perhaps consistent with this interpretation since the entropy does not 
change in these equilibrium models, and perhaps consequently supports the idea that P 2 represents a "gravitational 
entropy" . As we shall show W 2 and P 2 will be constant along homothetic trajectories, since all the dimensionless 
quantities remain constant along timelike homothetic trajectories. 
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A. The homothetic vector field 

The homothetic vector field (HVF) is calculated from equation 

£-Hogij = 2gij, (12) 

(see for example [13]- (4(| and [13]). Algebra brings us to obtain the following HVF: 

HO={t + to) dt + (l - (i + to) ~) xd x + (l - (t + t Q ) ^ yd y + (l - (t + t ) zd z , (13) 

with the following constrains for the scale factors: 

a(t) = a Q (t + t ) ai , b(t) = b Q (t + t ) a2 , d(t)=d Q (t + t ) a3 , (14) 
where 01,02,03 £ K, and the following restrictions for the constants 01,02,03 (obtained from the Eq. (|12p ) 

oi = 1, o 2 = a 3 . (15) 

As is observed we have been able to obtain a non-singular solution for the scale factors. Therefore the resulting 
homothetic vector field is: HO = (t + to) dt + (1 — 02) yd y + (1 — 02) zd z . 

Since we already know how the scale factors behave, then we may calculate all the curvature invariants as well as 
all the kinematical quantities. 

„ l + 2o 2 l-«i2 2 V&{a 2 -lf , , 

H = 77 - r , q = 2- — , a = — - T , (16) 

3 (t + t ) 2o 2 + 1 3 (t + £ )^ 

finding that 

= Lg ~ ^ = CO nst., (17) 
(l + 2o 2 ) 2 

where, as we can see, A £ (0, 1) , Va2 £ (0, 1) . Although the quantity A is constant, this quantity may take very small 
values, in fact A may runs to zero. It is also observed that the model never inflates since q £ (0, 2) , Va2 £ (0, 1) . 
Concerning the curvature invariants we find that 

Io = (6o 2 -2m 2 ) (t + t y 2 , 

h = 4 (3 (4 + ™ 4 ) - 4o^ + 2o^ (l - to 2 ) + 4m 2 (o 2 - 1)) (t + t y 4 , 
I 2 = 4 (2o 2 - 2o2 + 3fl2 - 2o 2 m 2 + to 4 ) (t + t y 4 , 
h = 16to 2 (to 2 + 6o 2 - 3 (1 + a%)) (t + i )~ 4 /3, 

h = 2to 4 (t + t Q y 4 /5, 

I 5 =2m 2 (a 2 -l) 2 (t + t y\ (18) 



and 



2 

TO 



(3 (1 + o|) + to 2 - 6o 2 ) 



W 2 = il = consi (19) 

9(l + 2o 2 ) 4 

, 4to 2 (to 2 + 6o 2 -3(l + o?)) , 

p 2 — i Z v ill = con st 20) 

3(2a 2 -2a3 + 3o 4 + (TO 2 -2a 2 )TO 2 ) 1 ' 

As above, although W 2 — const <C 1, dimensionless quantity, it may take very small values, for example, if 02 — > 1 and 
to^ (W 2 -> m 4 /3 6 ). 
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III. THE CLASSICAL MODEL 



We shall take into account the Einstein's field equations (FE) written in the following form: 

Rij - \Rgij = SvrGT™ - K 9ij , (21) 
where, A is the cosmological constant and T™, is the energy-momentum tensor defined by 

T ij = (Pm + Pm) UiUj + Pm9ij, (22) 

and where the 4— velocity is defined by: = (1,0,0,0), p m is the energy density and p m is the pressure. They are 
related by the equation: p m = ujp m , with lu e (—1,1]. In this section we study three models; vacuum solutions, a 
perfect fluid model and a model with a perfect fluid and where the constants G and A are time varying function. 



A. Vacuum solution 



In this case we have found only one solution a 2 = m = 0. Therefore, the metric Eq. ((BJ) collapses to this one: 

ds 2 = -dt 2 + (t + t Q f dx 2 + dy 2 + dz 2 . (23) 

We may compare this solution with the one obtained in the paper [24j where we were able to obtain a new solution 
belonging to Bianchi VI type. In this case, this solution does not belong to Bianchi VIq type, so we may say that 
the metric Eq. ([6j is more restrictive than the employed one in [24j . This solution is known as the Taub form of flat 
space-time (|25j chap. 9). 



B. Perfect Fluid 



For this model we obtain the following results: 



— -e 0,1, m =~J—± J -\ it..-. 

2 (w + l) V ' 2 V (tu + lf V 2 



(24) 



where u> € (— 4, l) , while the scale factors and the energy density behave as 

o(t) = a (t + t ) , b(t) = b (t + t ) a2 , d = do (t + i ) Q2 , p = po (t + i )" 7 , (25) 
where 7 = (1 + w) (1 + 2a 2 ) , po = g^j, ^4 = 2a 2 + a 2 — m 2 . With regard to the deceleration parameter 

1 (3cj -(- l) 2 

q = - (1 + 3w) > ^ = 1 — j = const. G (0, 1) , (26) 

Vo; e (— |, l) , so this model does not inflate, and A — >0 only when u> — > — |, while the Weyl parameter and the 
gravitational entropy behaves as 

W 2 = _ (3" + !)> -!) (*>+!) = cons< e (0j . 012) , (27) 
576 

P2 = 1 1 = const, e (-00, a], a^0 + , (28) 

3 (w — 1) (3uj z + 1) 

where as it is observed, W 2 < 0.01, it takes a very small values, W 2 <C 1, and it runs to zero if u — > 1 and ui — > —1/3. 
Notice that our solution is only valid if uj g (— |, 1J . Nevertheless P 2 has a very pathological behaviour. For example, 
P 2 -^a = 2.5000 x lO" 11 as uj -> -1/3, P 2 = when w = -|, and w = -| but P 2 -> -00 when w 1. 
Therefore, the metric collapses to this one: 

ds 2 = -dt 2 + (t + t ) 2 dx 2 + (t + t Q ) 2a2 (cosh 2mxdy 2 - 2 sinh 2ma;dydz + cosh 2mxdz 2 ) . (29) 
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Note that in [2J] we obtain two solutions, while with the metric Eq. ((6|) we are only able to obtain one solution which 
coincides with the one obtained in (24|. Therefore, this solution is valid when oj £ (— i, l) and m £ (0, -|] . It does 
not accelerate since q > 0. Nevertheless, we may say that the solution isotropizes since A — >0 when u> —> — | (then 
a-2 — > 1 = ai and m — > 0) and W 2 <C 1- The behaviour of P 2 shows us, that maybe, it is not a good definition for the 
gravitational entropy (at least in the framework of self-similar solutions) as we have already discussed. 

For ui = — |, and w = lwe get that m = 0, so the solution does not belong to Bianchi VIo type, furthermore if 
lu = l,then a 2 — i.e. we obtain the vacuum solution. 



C. Time varying constants model 

In this framework the FE are the following ones: 
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— SirGuip, 
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0, 

SnGujpi 
0, 



Ac 2 



Ac 2 , 



Ac 2 



A' = -8vrG'p m , 

where we have taken into account the condition divT — 0. 

Now, we shall take into account the obtained SS restrictions for the scale factors given by Eq. (fT5j). From Eq. 
we get 



(30) 

(31) 

(32) 

(33) 

(34) 

(35) 
(36) 



P = Po(t + t Q ) 7 , 

where 7 = (oj + 1) h and h = 1 + 2a2, since a\ — 1 and ai = 03. 
From Eq. (|30l) we obtain: 



A 



A (t + t y 2 - 8ttGpo (t + t y (uJ+1)h 



where A = 2a2 + a\ — m 2 . Now, taking into account Eq. (f31)| and Eq. (|3^|) . algebra brings us to obtain 



g = g (t + t y- 2 . 

while the cosmological "constant" behaves as: 



Gn — 



A 



Airpo (w + 1) ' 



In this case we have found the following solution 

a 2± = i (l ± - 4m 2 ) , 



Vm € 



1 1 

2' 2 



\{0} 



(37) 



(38) 



(39) 



(40) 



(41) 



0-2+ £ [1/2, 1), d2- € [0, 1/2), hence ft.± = 1 + 2a2± and therefore we obtain h + £ [2, 3) and /i_ € [1, 2) 

o(t) = a (* + t ) , b(t) = b Q (t + t Q ) a2± , d = d Q (t + to) a2± , 
P = po (t + t r J± , G = Go± (t + i ) 7± ~ 2 , A = A ± (t + t y 2 , 



(42) 
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with 7-fc = (lu + 1) h±. Notice that this solution is valid G (—1, 1]. In this way the metric collapses to Eq. (|29l) . 
The behaviour of the "constants" is the following one: 

decreasing if (lj + 1) h± < 2 [ < if (us + 1) h± < 2 

G«{ constant if (u) + 1) h± = 2 , A ± « < =0 i/ (w + 1) ft± = 2 , (43) 

growing if (w + 1) /i± > 2 I >0 i/ (w + l)/i±>2 



where /i± = l + 2a2± = 2±v / l — 4m 2 . Note that in 24| we obtained another solution. With regard to the deceleration 



parameter (for simplicity we have performed all these calculations with h+ i.e. with a2+) 

\W, (44) 



q + = 1 >0, q e (o, ] - ) v,n G 



1 1 

2' 2 



1 (v 7 ! - 4m 2 - l) 2 



A= - — % = const € (0, 0.06) -> 0, (45) 

4 (VI - 4m 2 + 2) 

, m 2 (4m 2 - 3 + 3Vl - 4m 2 ) 

W 2 = ; ; '- = const G (0, 0.0016 , (46 

18 (2 + VI -4m 2 ) V ; V ' 

„ 4m 2 (8m 2 - 3 + 3 VI - 4m 2 ) , x 

p 2 = v I i = const (47) 

3 ((6m 2 - 3) VI - 4m 2 + 12m 2 - 8m 4 - 3) 

where P 2 G (-co, 0.01) , Vm G [-|, ±] \ {0} . Therefore this solution is valid Vw G (-1, 1] and Vm G (-|, |) \ {0} . 
The model does not accelerate but isotropizes since W 2 —> as well as A <Cl. With regard to the quantity P 2 it is 
observed that P 2 — > 0, Vm G (—5, 5) \ {0} and only runs to minus infinity when m — > ±5. 



IV. SCALAR FIELD MODEL 



In this section we are going to study several scalar models. In the first place we study which kinds of potentials 
are compatible with the self-similar solution. For this purpose we study through the Lie group method the resulting 
Klein-Gordon equation. Once we have deduced the potential compatible with the self-similar solution we study if this 
kind of potential brings us to obtain self-similar solution. We answer in this case is no, we only obtain power law 
solutions but this fact does not mean that they are self-similar solution. Therefore, after this brief analysis on the 
potential, we start by studying a simple scalar model. In second place we study a non-interacting scalar and matter 
model. We leave for a forthcoming paper the study of the very interesting case of interacting scalar and matter models 
(see for example [HI and [1^]). In the third class of studied models we introduce the hypothesis of a G— var, i.e. we 
study a scalar model where G = G(t), is a function of time. In this case, in a phenomenological way, we outline a 
modified Klein-Gordon equation in order to take into account the possible variation on time of the function G(t). We 
go next to study the kind of potential compatible with a self-similar solution and G(t). To end this section, we study 
a model with scalar and matter fields and G— varying. 

The stress-energy tensor may be written in the following form: 

Z# = (p* + P+)ikUj+P*gii, (48) 
where the energy density and the pressure of the fluid due a scalar held are given by 

P* = l^ + V(<l>), P* = ^-V^), (49) 
while the conservation equation reads (Klein-Gordon equation) 

(/)" + <i)'H+ —V = 0. (50) 
a<p 

We need to study the class of potential compatible with the SS solution, for this reason we study by using the Lie 
group method the KG equation, (for an introduction to the Lie group method see for example [5(|-[53| and [54[ for a 
concrete application in cosmological contexts). In particular, we seek the forms of V(tfi) for which our held equations 
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admit symmetries and therefore they are integrable. In this case we already know that the Hubble function behaves 
as: H = h(t + <o) _1 i h G M + , so the KG equation reads 

0" + h^'(t + t y 1 + ^ = o. (51) 

a<p 

Theorem 1 The only possible form for the potential V (</>) for a spacetime admitting a HFV, HO, is V(</>) = 
Vb exp (Kip) and therefore cf> = alnt. 

Proof. The application of the Lie group method brings us to outline the following system of PDEs 

t^t 2 = 0, (52) 
2htfy + t 2 ^ - 2t 2 i t4 , = 0, (53) 

3t 2 ^ — + htit -H + 2t 2 rj tlP -t 2 £ tt =0, (54) 
d(p 

2 d 2 V 2 „ 2/ dV 2 dV n 
d<p' ! a<p d<p 

If we impose the symmetry £ = a (t + to) , r] = 5, then its invariant solution is: = — In — (t + to) , then, from Eq. 
([55]) . we obtain the next restriction for the potential V 

d 2 V „ dV . „ „ / a 



^——75- + 2a— = => V = / 9exp -2 T </» + k, a, /3, 8, k € M. (56) 
dcp* d<p \ J 

Therefore we have found, redefining the numerical constants, that the only solution compatible with the FE is 

cf) = ±y/aln{t + t ), V = (3 exp (^-^=^j . (57) 

as it is required. ■ 

Note that in this case it is possible to find more symmetries, but the solution generated by them are not compatible 
with the FE. For example, if we impose the symmetry, £ = at, r\ — 5<f>, then Eq. (155 j) yields 



+ ( 2a ™ <5) ?T = ° => ^ = «if |M) +«2, (58) 
d<p z d<p 

which is the potential proposed by Peebles and Ratra, V « 4>~ a 0, but this solution it is not compatible with the 
FE with a SS solution. Nevertheless we shall use this potential in the G— varying scenario. In the appendix we give 
an alternative derivation of all these results by using the matter collineation approach following a previous paper (see 
.54]). 

Models with a self-interaction potential with an exponential dependence on the scalar field of the form V = 
f3 exp (=F2</>) , have been the subject of much interest and arise naturally from theories of gravity such as scalar-tensor 
theories or string theories [55| . Recently, it has been argued that a scalar field with an exponential potential is a strong 
candidate for dark matter in spiral galaxies [56| and is consistent with observations of current accelerated expansion 
of the universe [53] . 

In the inverse way we may state the following theorem. 

Theorem 2 For a scalar model if the potential is of the form V = f3 exp (t2</>) , then the scale factors must follow a 
power law solution i.e. H = ht^ 1 . 

Proof. As above we perform the proof by using the Lie group method. In this case we must study the following ODE 

(f>" + <P'H+-^V = 0, 

where V = /3exp (=F2</>) , then we shall study the different forms for the function H(t) in order to get and integrable 
ODE. 
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We have the next system of PDEs 



St* = 0, (59) 

+ m - = 0, (60) 

-Ge- 2 % + H& - H'Z + 2 m<j> - £ it - 0, (61) 

4e- 2 *77 + ijtt - 4e-*+£ t + 2e~ 2 % + H m = 0, (62) 



As we can easily see, the symmetry £ = t, r\ = 1, brings us to get = lnt, as invariant solution, and from Eq. (joTjl we 
obtain the result i.e. if = ft. G R. ■ 

Obviously this result does not mean that the solution must be self-similar, it only means that the scale factor must 
follow a power-law solution i.e. they are of the form a,(t) = aot aj , with dj € M + . 



Scalar model 



We write the FE in the following form: 



ab ' a 7 + 7I-( 2+ - + -)- =87rG ^ (63) 



a' b' a' d' d' b' 
b ad d b 
b" d" d! b' 
t 

a" b" a' b' 

T + T + 76~( 2 + ^^)^ = ^ Gp ^ (65) 



b 'T + dT + ( 2 + ^ + ^'~ = _8 ^' (64) 



6 d a 6 a d 



b 2 


d 2 


\ TO 2 




b 2 


J 4a 2 


6 2 


d 2 


\ TO 2 


f d* + 


b 2 


J 4a 2 


3d 2 


lr 


\ TO 2 




d 2 


J 4a 2 


( 62 






V d 2 a 2 


l 


i 2 a 2 ) 



~J + ~ + ~a~~d~ ( + ~d 2 ~~b 2 )Aa^ 2 



m Mi-I^)=0, (66) 
8ttG P ^ (67) 



b' 2 H+^-V = Q. (68) 



By assuming the potential given by Eq. (|57| it is possible to find the next set of solutions 

\{0}, 



«2± = ^ (l ± \/l - 4m 2 ) , Vm G 



1 1 

2' 2 

a(t) = oo(t + t ), b(t) = b (t + t ) a2± , d = d (t + t ) a2± , (69) 

and 

a± = 1± \/i -4m 2 , /3± = i (l ± v 7 ! - 4m 2 ) 2 , (70) 
= ±V5±Mt + to), ^ = /3±exp (t-J= <A ■ (71) 



As it is observed, we have obtained the same behavior for the scale factor as the one obtained in the case of a perfect 
fluid with time- varying constants model. For this reason, as we already know, we get: q > 0, Vm G [— i, i] \ {0} , 
.4. = const G (0, 0.06) — ► 0, while the Weyl parameter and the gravitational entropy behaves as W 2 = const g 
(0,0.0016) <C 1, and P 2 — const e (— oo,0.01), Vm e [— |, |] \ {0} . Therefore the model does not accelerate but 
isotropizes since the quantities (A and W 2 ) instead of being constant, they take values very close to zero. With regard 
to the quantity P 2 it is observed that P 2 — > 0, (it takes values very close to zero) Vm g ( — 5i f ) \ {0} an d it only runs 
to minus infinity when m — > ±«. 



B. Non-interacting scalar and matter fields 

The stress-energy tensor may be written in the following form: T = T m + T^, where the energy density and the 
pressure of the fluid due a scalar field are given by Eq. (|48l) . This describe a non-interacting dark matter and dark 
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energy cosmological model (we assume that the baryon component can be ignored). Since the nature of both dark 
energy and dark matter is still unknown, there is no physical argument to exclude the possible non-interaction between 
them. 



a 1 V a 1 d! 
a b ad 



We write the FE in the following form: 

u # d' b' 
~d~b 

w_ <f_ d/_ v 
T + ~d + IT 

a" b" a 1 V 

— + -r + —r- 

a b a b 



~d? 



d 2 \ m 2 

2 



b^__dr_ d_V__d_di 
b d a b ad 



b 2 d 2 \ mr 
2 H 1 

d 2 b 2 4a 2 



3(P \P_\rrf_ 
l^^d 2 )^ 2 
b 2 d 2 



d 



-r + — + 7 - 2 



a!_d!_ 
a d 



d 2 c 
d 2 



b 2 c 
b 2 } 4a 2 



8nG (p m + 
—8ttG {ujpn 
—8irG (ojp„ 
0, 

—8ttG (ojpn 



and the conservation equations now read 



and 



p' m + (w + 1) Pm H = 0, 



'H+—V = 0. 

dcp 



where H = h (t + to) 1 , and h = 1 + 2<22. 
In this case we have found the next solutions 



«2 



l-U) 

2w + 2 



Ckd2 



e(o,i) 
( 1 -<■ 



\2uj 



1 



Po 



1+UJ 

■ uj — a (1 + u>) 



1 y/- (3^ + 1) [u - 1) 
m = t: : € 0, - 

2 1 +uj \ 2 



(1 



>0, 



(72) 
(73) 
(74) 
(75) 
(76) 

(77) 
(78) 



(79) 



where the constant a must verify the condition: a > > 0, £ (— |, l) . Therefore we have the following 

behaviour for the main quantities 



ax =a (t + t Q ) , b = b Q (t + t ) a2 , d = d (t + t ) c 



(80) 
(81) 



iy/a In (t + to) 



V = $ exp ( =F 



(82) 



Since the scale factor behaves as in the perfect fluid solution (see above) then the deceleration parameter (as we 

already know) behaves as: q = | (1 + 3w) > 0, A = — const. £ (0,1), Vw e (— -|,l). Therefore, with 

the above restrictions on the oj— parameter our model does not inflate, q > 0. While the Weyl parameter and the 
gravitational entropy behave as: W 2 — const £ (0,0.012), and P 2 = const, £ (-co, a], with a — > + (see the above 
discussion about these quantities). 



C. G— varying 



We would like to study how the gravitational varies constant when we are considering only a scalar field. For this 
purpose, in analogy with the perfect fluid case (see [Uj]) and in a phenomenological way, by using the Bianchi identity 
div (8nG(t)Tij) — 0, we propose the following conservation equation 



Gp' + G (p + p) H = -G'p 



(83) 
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which is the modified KG equation. 
For this model the FE read 



a'b' a'd' d' V ( b 2 d 2 



ab' ad' db V 2 + ^ + ^J^ = 8 ^ W ' (84) 
b" d" d'V / b 2 d 2 \m 2 

T + l[ + l-b + [ 2 + lp + lp)^ 2 = - 87rG(<) ^' (85) 
a" b" a 1 b' ( U 2 b 2 \ m 2 
- + T + -a-b-{ 2 + l?-d 2 )^ 2 - = -^ G ^ 
b" d" a' b' a! d! 2 ( b 2 



b dab ad' \d 2 a 2 b 2 a 2 ) ' ^ ^ 

d" a" a' d! ( 3b 2 d 2 \ m 2 

In order to solve the FE we need to solve Eq. (jg3"j) . To that end, we shall study it through the LG method. Eq. 
(|83| could be rewritten in the following form 

^V + ^-V 2 + ^' + P^=0, (90) 

where H = ht^ 1 . For simplicity, and without lost of generality, we consider H = ht^ 1 instead of its non-singular form. 
As above, we are seeking the forms of V (cj>) and G(t) for which our field equations admit symmetries and therefore 
they are integrable 

As above, in order to study the possible solutions to Eq. (|90|) we apply the LG method, where the standard 
procedure brings us to get the following system of PDEs 

Z<p<p = 0, (91) 
V4>4> ~ + Ufa' 1 + = 0» ( 92 ) 

2Vt* +(\{^-~ (§) j - ht ^ £ + 3 ^ + (to- 1 + 6 - ft* = 0, (93) 



??tt + V^r) + 4-^- + (At" 1 + —J ry* + 2V*& " V^Tfe = 0, (94) 

Gil / /nr/ \ 2\ /^f/ /nr/ /nr/ 

/ \ \ L» 



G \G) i n + ~G V4,V + 3 ~G nt ~ 2 ~G Vr]4 '^ ' (95) 

= 0, (96) 



where, V<f, = 4j. The following symmetry 



e = — , »? = ==► = (97) 

a 

then we obtain the following restrictions, from Eqs. (I93ti95[) . From Eq. (|93|) we get 

G " = 7T~T ^ G = «i« 9 , (98) 

while from Eq. (jMf we obtain 

VW^-f- + lV*=° =>■ ^ = ^ 2 (^ +1 ), (99) 



12 



in this way V = K 2 (t + t y 2(a+1} . So we have found that the main quantities behave as follows 

<t> = (t + h)- a , V = n 2 ^ +1 ) =K 2 (t + i r 2(a+1) ,G=Ki(t + t ) 9 , (100) 

such that g — 2 (a + 1) = —2, and therefore g = 2a. Note that we may redefine the constants in order to get V ~ <\T a . 

Theorem 3 The only compatible form for the potential V ((/>) with the FE for a spacetime admitting a HFV, HO, 
where G — G(t), is V(<f>) — Vo<f)~ a and therefore 4> = (t + toY , and G = K\ (t + to) 9 , with a,j3,g£ M. 



Taking into account all these results we have found the next solution 



a 2 = - (l ± y/l- 4m 2 j , Vm e 



1 1 

2' 2 



\{0}, (101) 



1 - 2m 



2\ 



a = a> 0, G = const. > 0, (3 = ± — — '- + y/ (I - 4m 2 ), 

Go 

so the behaviour of the main quantities is the following one 

oi = ao (t + t ) , b = b (t + tof 2 , d = d Q (t + tof 2 , (102) 

$ = $o{t + toT a , V = P(t + to)- 2{a+1) , G = Go(t + t ) 2a . (103) 

Notice that this is the same solution for the scale factor than in the above models. Therefore, with the above 
restrictions on the m— parameter our model does not inflate q > 0. i.e. the model does not accelerate but isotropizes 
since W 2 — > as well as A. With regard to the quantity P 2 it is observed that P 2 — > 0, Vm e (—5, |) \ {0} and it 
only runs to minus infinity when m ±5. With regard to the gravitational constant G, we have obtained that it is 
an increasing time function since a > 0. 

D. G variable with matter and a scalar field 

We start by rewriting the stress-energy tensor in the following way, T 1 ^ = + , where T y = (p + p)u % u? J rpg L ^ , 
and p = p rn + ptf, and p = p m + p^, and taking into account the Bianchi identity 

div (8nG(t)Tij) = <^ Gp' + G (p + p) H = -G'p, (104) 



i.e. 



/ \ Q) 

P'm + (Pm+Pm) H + <j)' \U4> + —j =-— ( Pm + : 



(105) 



We may study Eq. (|105[) in several ways. One of them, maybe the simplest one, may be spliting it into 



Qt 

P'm + (Prn+Pm) H = ~—p m , (106) 

<io7 » 

Notice that this approach is similar to a scenario describing an interacting scalar and matter fields. 
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There the FE for this model are: 



a'b' a' d! d! V / b 2 .. , ,,. 
b" d" d'b' ( b 2 



b" a'b' / 3d 2 



a 



- + T + - I - 2 + -- - - = -SnG(t) (u Pm+P ,) , (110) 



b" d" a' b' a'd' 2 ( b 2 d 2 



b d • ab ad < m {d^-^J = °' (1U) 
d" a" a'd' / 3b 2 d 2 \ m 2 . 

G" 

Pm + (Pm+Pm)-ff = --gPm, (113) 

=4^G* ,2+,,( 4 <ii4) 

With a potential given by Eq. (|100[) we have found the next solution 



a-2 



- (l ± \/l - 4m 2 ) Vm G 



2 



1 1 
2' 2 



\ {0} , a = a > 0, 



= 



(l ± ^/(l - 4m 2 )) + G a (1 - w) + 1 ± - 4m 2 ) - 4m 2 (1 + w) 

2G (w + 1) : 



thus 



1 ± ^(1 - 4m 2 ) - G a 2 „ , s 

Po = n ! Til — ' G ° = const > °' 115 

G (w + 1) 



Pm = Po(* + *o) (( " +1),1+2q) , p m =U>p m , (H6) 

oi = ao(t + to), 6 = &o (t + ioP , d = d (£ + to) a2 , (117) 
4 = 4>o (t + t y a , V = (3 (t + t )" 2(Q+1) , G = G (i + t ) 2a . (118) 



Therefore, as in the last model, with the above restrictions on the m— parameter our model does not accelerate but 
isotropizes since W 2 — > as well as A. With regard to the quantity P 2 it is observed that P 2 — > 0, Vm € (— |, 5) \ {0} 
and only runs to minus infinity when m — > ±4. G is an increasing time function as in the above model. 



V. SCALAR-TENSOR MODEL 

We consider the following field equations for the BD model (59j . 

Rit - \ 9ij R = jT™ + A ^ + -| Lrfj - + ~ (ct> ;ij - g^Ucf) , (119) 

1 I d ^ OJ (0) ^1 { n , n d 



^ + 2^W n {^) + 2^W){ R + 2 dt ^ M ) = °- < 120 > 

The arbitrary functions uj (</>) and A (</>) distinguish the different scalar-tensor theories of gravitation. A (0) is a 
potential function and plays the role of a cosmological constant, and to ((/>) is the coupling function of the particular 
theory. is the matter stress-energy tensor. 
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The last equation can be substituted by 

D * + jrhm - * A <*>) ~ wrhrn ( 8 * T - • (121) 

where T = T\ is the trace of the stress-energy tensor, where we have assumed <j> — <j>(t), and the derivatives respect t 
are denoted by a comma. Furthermore it is verified the following relationship: divT = 0, i.e. 

p' + {p +p)H = 0. (122) 

In what follows we shall assume to ((f)) — const, A — A (ft. The corresponding field equations with a perfect fluid for 
the matter content in the homogeneous line element (Bianchi VIo model) will be calculated. Thus the field equations 
are 

d 2 \ m 2 8tt TT ft lu fft^ 2 



a' b' a' d! 


d! b' 




b 2 


a b ad 


' ~d~b ~ 


(- 


' d 2 


b" d" 


d! b' 




b 2 


~b + ~d~ ] 


' d b + 




' d 2 


a" b" 
a b 


a'b' 
a b 1 




3d 2 
b 2 



b 2 )^-^- H -6 + 2{V +AW ' (123) 



2 



-P- I l\^ + ^)--[ I i) -^+A(0), (124) 



n 2 _ 8tt ft fd' , b'\ to (<f)' 
a 2 

b 2 \m 2 8tt ft fa' b'\ oofftV ft' A , . , 2 ft ( d! b 

(125) 

b~ d + a{b~ d) + a 2 [d 2 ~ b 2 ) ~ cf>{b d)' 1 ' 

d" a" d d! ( 3b 2 d 2 \ m 2 8ir ft f a 1 d'\ cofftV ft' WN , , ,2 ft ( d' b 
^ + - + ^l-\ 2+ ^-Y 2 )^ 2= -^ P '^[^ + l)-2\^) - T +A(0)-cosh(mx) ^^-- b 

(127) 

(3 + 2. (ft) + h€) - 2 (A - ^) = * (p - 3p) , (128) 

p' + (p+p)tf = 0. (129) 

Since we are only interested in finding self-similar solutions then, if we take into account our previous results, i.e. 
b = d, the FE reads 



4^ 




)'- 


2 

m 














)'♦ 


2 

m 
'a 7 


a" 6" 


a' 


&' 


m 2 


— + T H 
a o 


a 


& 


a 2 


d" a" 
a a 


a' 


d' 


2 

m 


a 


d 


"a 2 



8tt rr 0' w /</>' N 2 

_8tt w ' ' 2 

~ (f) P ~ (f> b~ 2 

8tt ft fa' b'\ co fft^ 2 



(130) 
(131) 



*>-j{a + T)-2{j) -T +A(0) ' (132) 
8tt 0' fa' b'\ to fft^ 2 



-P-j{a + J -* [ ^ ] -^ +A(0) ' (133) 



and the conservation equations 



(3 + 2oo(ft))(^-+H^-) -2(A-ft^-) =^(p- 3 p), (134) 



p' + (p + p)if = 0, P = Pat- a , (135) 

where i? = /i (i + to)~ , with ft- = (1 + 2a 2 ) , and a = h (1 + 7), we are taking into account the equation of state 

p = 7p, 7 e (-1, 1]. 

In order to solve the resulting FE we need to integrate 

>(£+^) ->(*-*£) -7^- 
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We may study this equation through the LG method, i.e. we study the kind of functions A (0) such that this 
equation is integrable in a closed form. We start by rewriting it in an appropriate way 

4>" + hr 1 ^ - B ( A - <f>^j) <j> - Cr a = 0, (137) 



where 



d(j>J 



and we shall deduce that a = 2 — n. 

We need to solve the following system of PDEs 

< 2 £a<A = 0, (139) 
2ht- 1 ^ + 77^ ~ 2£ 0t = 0, (140) 
-3 {B<j> (A - 0A ) + Cr a ) ^ + ht~ 2 {ti t - + 2% - Ztt = 0, (141) 

B V (tfA^ - (A - 0A )) -250 (A - 0A ) £ t + -^ f a| - 2&J + (s^ (A - + C^) v^-m+Vu = 0. (142) 
For example, if we impose the symmetry 

£ = t, r] = n<t>, (143) 
brings us to obtain the following restriction on A (</>) . From Eq. (|142[) we get 

B<t> (n^A^ - 2 (A - 0A^)) + Ci~ a (a - 2 + n) = 0, (144) 
and therefore, we find as result, that: 

n = 2-a, (145) 

and 

n</> 2 A 00 - 2 (A - 0A ) = A = <f>~ 2 / n , (146) 

is a solution. In fact the most general solution is A = C\4> + C2<j>~ 2 t n ■ This result is valid for all the self-similar 
Bianchi models. Furthermore, the symmetry Eq. (|143l) brings us to obtain a particular solution of Eq. (I137[) which 
is given by, <f> — 4>ot n , with (f>o £ K, with n = 2 — a, and the following constrains on the numerical constants must be 
verified: n(n - 1) + hn - 2(1 + %)B - C = 0, a = h (1 + 7) = (1 + 2a 2 ) (1 + 7) . We would like to emphasize that 
other solutions could be obtained with this procedure by imposing other symmetries. 
We find the following solution: 

<Po = 4>o 4>o = 1, A = 0, 

( 7 -l)M 7 -l)-l) 
G2 = - 2c( 7 2 -l)+ 7 -3 ' a2= °' 7 = 1 ' 

2 (a; ((37+I) (7-1)) -2) 
9= 4c( 7 -l) + 3 7 -5 ' 9 = °' 7 ^ ±! 

(3 + 2c) 2 ( 7 -l) 3 

= (7 - 1) V-2 (3 + 2c) (a, ((3 7 + 1) (7 ~ 1)) - 2) 
2 (2w ( 7 2 - 1) + 7 - 3) 
m = 0, 7 =1AA ± , 

- -(37-D(7-l) „ = 0i _ 7 _ 1A 1 



2w( 7 2 - l) + 7 -3' ' ' 3' 
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where A± = 5 -. 

We get for the BD parameter oj that according to solar system experiments is oj « 500. (see [60]). A better 
estimation of this parameter should be obtained from measure of other cosmological parameters in order to constrain 
oj more strongly than by means of solar system experiments (see [ll|). However, theories of the very early Universe 
such as string theory, are better described in the context of JBD, which shows that oj can take negative values (see 
for example [Uj]). A recent value for oj is oj « 3300 (l6j . 

If we fix oj — 3300, then we get the following results 

A+ = 1.0002, = -0.33348, m > 0, V 7 G (A-, 1], 

( <0 V 7 <i_ 
o 2 >0, V 7 G (-1,1], q= { =0 7 = A_ , 

[ > V 7 > A_ 

( > V 7 < 1/3 

Po>0, V 7 G(-1,1], n= =0 7 =1/3A1 . (148) 

[<0 V 7 €(|,l) 

Therefore this solution, with oj — 3300, is only valid V 7 G 1]. Note that if 7 < A_ then to is not defined. This 
means that 

a = a (t + t ) , b = b (t + t ) a2 , d=d (t + t ) a ' 2 , H = h(t + to)' 1 , 

P = Po{t + t Q r a , (j, = (t> (t + to) n , A = 0, (149) 

and therefore the scale factors are increasing time functions, the energy density is a positive time decreasing function. 
The solution does not inflate since q > V 7 £ (A-, 1]. <fi is a positive growing time function if 7 G (A-, 1/3), it 
is constant if 7 = 1/3 and it behaves as a decreasing time function if 7 £ (1/3,1), if 7 = 1 then it behaves as 
a constant. This means that G is a decreasing time function if 7 G (A_,l/3), it behaves as a true constant if 
7 = 1/3 A 1 and if 7 G (1/3, 1) then it is a growing time function. The cosmological "constant", A vanishes. Notice 
that a = (l^r^l^Ka+ii 

21^(7^ — 1)+7— 3 

In order to check if the solution isotropize we compute the quantities A and W 2 , obtaining 
A (w( 7 2 - 2 7 -l)-2) 2 

.4 = -^-^ '- V = const G (0, 1 , (150 

(4w( 7 -l) + 3 7 -5) 2 



w 2 = - (7 J?. 2 ( " (7 ' V r ( 8w ((27 + 1} (7 " 1)} " 372 + 67 " 15 ) ' (151) 

36 (4w ( 7 — 1) + 3 7 — 5) 



W 2 = const G (0,0.01), V 7 G (A-,l\ and oj = 3300. -4(A_) = W 2 (A_) = 0. With regard to the gravita- 
tional entropy we have obtained the following behaviour: P 2 (A_) = 0, if 7 G -0.2000969530) then P 2 > 0, 
P 2 (-0.2000969530) = and it runs to -00 V 7 G (-0.2000969530, 1], so once again, we have checked that this quantity 
is not a good definition for gravitational entropy. 
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A. The particular case 7 = 1/3. 
In this case T = and therefore we find the next solution 



9o 
A 

0,2 

q 

Po 

a = 



90, 
1 
6 

3 V 



1, 



(3 + 2u) (Am + Vtj 
( 3m + Vs) , a 2 = 



A o = 



m — — - 



m = — - 



V3 
3 ' 



-2m 



2m + V3 
1 

~16^ 
4^ 



q = 



32a; m 



> Vra < 
< Vm > ' 

+ 9 + 24toV3 + 44m 2 ) , 



(2m + Vtj 



Po = 



in = — 



n = — - 



(4m + Vtj 



V3 (4w ± V2w + 3 + 6) 
16u; + 22 ' 

n = <=^> m = — 



V3 
4 ' 



If we fix cj = 3300, then we get the following results 

Po > 0, Vine/, a > 0, 
A > 0, Vm e /, A = m Ao 



ft2 > 0, Vm G /, a 2 = m = 

g > 0, Vm e /, 

> Vm e (-0.43569, m Ao ) 
n = ( = m = mA 

< Vm e (m Ao , -0.43036) 



Vm G /, 
4 



—- = -0.43301, 



(152) 



(153) 



therefore this solution, with to = 3300, is only valid Vm e I, where I = (-0.43569, -0.43036). Note that m Ao G /, 



This means that 



A (mA ) = °i fl 2 ( m A ) = 9 (™a ) = 1) Po (mA ) = 1-4921 x 10 2 . 



a = a (t + t ) , b = b (t + i ) a2 , d = d Q (t + t ) a ' 2 , H = h(t + t ) 
p = Po (t + toy a , 9 = 90 (t + toT , A = A (t + t ) -2 , 



(154) 



(155) 



and therefore the scale factors are growing time functions, the energy density is a positive time decreasing function. 
The solution does not inflate since q > Vm G /. (f> is a positive growing time function if m G (—0.43569, mA ), it is 
constant if m = mA and it behaves as a decreasing time function if m G (m Ao , —0.43036). This means that G is a 
decreasing time function if m G (—0.43569, m\ ) , it behaves as a true constant if m = m Ao and if m G (—0.43569, m\ ) 
then it is a growing time function. The cosmological "constant" , A is a positive decreasing time function except in 
m = m Ao . 



VI. CONCLUSIONS 



In this paper we have studied some Bianchi types VIo (with an unusual metric) models under the self-similarity 
hypothesis. We have started by comparing our results with the "classical" perfect fluid solution already studied by 
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Collins, Wainwright and Hsu and other authors [24[. Furthermore, we have been able to improve the solutions since 
we have found a non-singular solution for the scale factors i.e. they behave as a(t) ~ (t + t ) ai . However, the metric 
employed in this paper is very restrictive, since it allows us to obtain less solutions than with the usual one [24l |. 
Nevertheless we have been able to obtain a new solution for the case of a perfect fluid with time-varying constants. 
This solution is not inflationary but it is very close to isotropizing since the quantities A and W 2 take values very 
close to zero. In fact, for an adequate selection of the parameters u and m they run to zero. This solution is valid for 
all uj G (—1, 1] and m G [—1/2, 1/2] . In this case we have been able to enlarge the range of validity for the equation 
of state and we have shown that if G behaves as a growing time function then A is a "positive" decreasing time 
function. In the same way, if G is decreasing then A behaves as a "negative" decreasing time function. With regard 
to the gravitational entropy, we have come to the conclusion that the quantity P 2 is not an acceptable candidate for 
gravitational entropy along the homothetic trajectories of any self-similar spacetime (in all the cases studied in this 
paper). 

In the second model we have studied a scalar field. We have started this section by calculating the potentials 
compatible with the self-similar solutions. Inversely, we have proved that for such scalar fields the scale factor must 
follow a power law solution. These theorems are very general and are valid for all Bianchi models. We have studied 
two cases. From the first one, with a scalar field alone, we have obtained a solution that is not inflationary but it could 
be considered to be very close to isotropize, since, as above, since the quantities A and W 2 take values very close to 
zero. In the second case, we have studied a non-interacting scalar and matter fields. The solution is not inflationary 
but isotropize as in the previous cases, and it is valid Vw G (—1/3, 1) and m G (0, 1/2]. In order to incorporate into 
this framework a variable G, we have proposed, in a phenomenological way, a modified Klein-Gordon equation. We 
have studied the kind of potential compatible with a self-similar solution and a variable G. Once we have deduced the 
potential and the scalar field then we study two cases, a scalar field with a G— var and a scalar field with a matter field. 
The conservation equation outlined in this case is quite similar to the one employed in the case of interacting scalar 
fields. The solutions obtained are similar since the scale factor is the same and therefore they are not inflationists 
and close to isotropize. In both cases, G behaves as a positive increasing time function. 

In the scalar-tensor model, for simplicity, we have chosen, lj (</>) = const, and A (0) playing the role of an effective 
cosmological constant. As we have shown, the resulting FE are quite difficult to study. Nevertheless, since we are 
only interested in studying self-similar solutions, we have been able to simplify the FE. We would like to stress that 
we have not needed to make any assumption in order to integrate them. By using the Lie group method we have 
obtained a possible form for the dynamical cosmological constant, A (0) = 4>~ 2 / n — t~ 2 , since 4> = (f>ot n . In the same 
way, we emphasize that this result is valid for all the self-similar Bianchi models. With this result we have obtained 
several solutions for the model. We have considered that the first of the obtained solutions is unphysical since po < 
if the coupling parameter lu is positive as the recent observations suggest. Nevertheless if we consider u < 0, as it 
is suggested by the string theories (62[, this solution has physical meaning. In the second solution, A = and it is 
only valid when 7 G I, where / = (A-, 1], where A_ = —0.33348 if oj = 3300 as recent experiments suggest [la ]. For 
such values </> behaves as a growing time function if 7 G {A—, 1/3) , it is constant if 7 = 1/3 A 1 and it is a positive 
decreasing time function in the interval 7 G (1/3, 1) . Therefore, G is decreasing, constant and growing in the same 
intervals. In the same way we have found that this solution does not inflate since q > 0, 7 G /, which is unusual. We 
may also say that this solution would be considered to be very close to an isotropy state since the Weyl parameter, 
W 2 -C 1. In fact this quantity takes values very close to zero 7 £ J. We have also studied the particular solution 
7 = 1/3. In this case the trace of the stress-energy tensor vanishes and therefore the conservation equation reduces to 
a very simple ODE. This solution is only valid for a very restrictive interval, m G I m = (mi, 1712) , where po > 0. We 
have found that in this case A behaves as a positive decreasing time function except when m = m\ a G /, for which, 
A (toa ) = 0. In the same way, we have found that ^ is a growing time function when m G ™a ) , it is constant 
if jn — m\ and it behaves as a decreasing time function if m G {m^ a , 7722 ) . G behaves in the inverse way in the same 
intervals. This solution does not inflate as the one above does. Finally we would like to stress the fact that, in this 
solution, the exponents, 02, n and a are irrational numbers. This could have implications for the integrability of the 
FE (see [gl, HI]). 

Appendix A: Matter collineation for the scalar model 

In this section we shall study the matter collineations for the scalar field following a method employed in [HJ]. We 
start by defining a generic vector field X = (X. t (t, x, y, z))? =1 G X(M). The energy-momentum tensor is defined by Eq. 
(1221) . The metric tensor g^ is defined by Eq. (|r?j). In recent years, much interest has been shown in the study of matter 
collineation (MCs) (see for example |65|-[73|). A vector field along which the Lie derivative of the energy-momentum 
tensor vanishes is called an MC, i.e. CxTij = 0,where X 1 is the symmetry or collineation vector. Also, assuming 
the Einstein field equations, a vector X 1 generates an MC if CxGij = 0. It is obvious that the symmetries of the 



19 



metric tensor (isometries) are also symmetries of the Einstein tensor Gij, but this is not necessarily the case for the 
symmetries of the Ricci tensor (Ricci collineations) which are not, in general, symmetries of the Einstein tensor. If 
X is a Killing vector (KV) (or a homothetic vector), then CxTij = 0, thus every isometry is also an MC but the 
converse is not true, in general. Notice that collineations can be proper (non-trivial) or improper (trivial). A proper 
MC is defined to be an MC which is not a KV, or a homothetic vector. Carot et al (see |66|) and Hall et al.(see |67[) 
have noticed some important general results about the Lie algebra of MCs. Let M be a spacetime manifold. Then, 
generically, any vector field X e X(M) which simultaneously satisfies LxTab = (<<=> LxG a b = 0) and Cx^Scd = is 
a homothetic vector field i.e. Cxg — 2g. 

The usual matter collineation equations read, CxTfj = 0, Tf- is given by Eq. (|48p. finding, in this case, that the 
obtained matter collineation (MC) is: 



X = Xidt + (X[ - XxH 2 ) yd y + (X[ - XiH 2 ) zd Zl 



(Al) 



which is a proper MC, and where, as it is observed, if X\ = (t + to), then it is regained the usual homothetic vector 
field (see Eq. (O) i.e. a improper MC. 

For this reason we may also check that the homothetic vector field also verifies the equation, CnoTij = 0, (note 
that it is verified CnoC^ cd — 0) that we may develop as follows: 

p't + 2p = 0, 



.9332/ 



.9342/ 



b' b" b' 2 

~b +t T~ V 

b' b" b' 2 

b +t T~ V 



aa + taa" — t (a 1 ) = 0, 
(d 1 d" d' 2 \ 

(d! d" d' 2 \ n 



322 ( tp' + 2p- 2pt- ) + tpg' 22 = 0, 



<7 33 ( tp' + 2p- 2pt- 

, , ;v d 1 

5 34 ( tp' + 2p - 2pt ( - + - 



544 tp' + 2p- 2pt 



d! 



As we can see, actually, the only ODEs that must by satisfied are: 

p't + 2p = 0, (p't + 2p) = 0, 

which are equivalent. Hence 

p't = -2p, L H p = p't = -2p 



0. 



tpg'33 + Pxd x g33 ( 1 - t— 
tpg' 34: + pxd x g 34 (l-t^) = 0. 
tpg'44 + pxd x gu ( 1 



(A2) 



(A3) 



(A4) 



i.e. 



that we may split into 



dV ^^ t +2 (V + ^)) 



t (4,"4>') + 0' 2 = 0, t^-4>' + 2V = o 



(A5) 



(A6) 



where 

t(f>" + (j)' = 

i.e. Ch4>' — 0. With regard to the second equation 

dV 



k + 2V = 



K In i, 



V = Ke~-' 



(A7) 



(A8) 
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i.e. C H V = -2V. Note that from Eq. ([A7[l. fit = k. 
We also may study the complete equation (IA5|) i.e. 



0" + ^ + 0'^ 1 +2^(^') _1 =0, (A9) 
d</> 

under the LG method. The standard procedure brings us to get the next system of PDE: 

t 2 ^ = 0, (A10) 

t 2 n H - 2t 2 £ t4> + 2*6,5 = 0, (All) 

2t 2 n t ^ - t 2 itt + t^ + St 2 ^^- - £ = 0, (A12) 

acp 

t 2 m + St^V + 2t 2 Zt^ - 2t2 ^% + tr lt + ^~^ = °i ( Al3 ) 

-2^ + 6i&7 + 2ir?^ - 4%V = 0, (A14) 
d<p 

tri t V = 0. (A15) 
The symmetry, £ = at, rj = 5, brings us to obtain the following restriction on the potential (from Eq. (|A13|) and 

dV d 2 V dV , kS 

2 ^ + w =0 > 2t/+ ^r°< (A16) 

and therefore we obtain as solution 

V = cxp (~2fi) (j) = hxt. (A17) 

Therefore we may state the following theorem. 

Theorem 4 The only possible form for the potential V {fi) for a spacetime admitting a HFV, HO is V{fi) — 
Vq exp {nfi) and therefore (j) — hit. 

Sometimes it is interesting to study the symmetries of the tensor T- £ Ti(M). In this case the matter collineation 
equations read ChoT^ = 0, iff p' = 0, and p' = 0, which is equivalent to 

(A18 , 

where as we can see this approach is related with the variational symmetries. 
In this case the solution of Eq. (|A18[) is the following one 

t = [ ± da + C 2 , (A19) 



y/-2V(a) + d 
The Lie group methods applied to Eq. (|A18|) gives 

= 0, (A20) 

VH> ~ = 0, (A21) 

2^-6* + 3^^=0, (A22) 

„ ,dV dV d 2 V „ , A „ oN 

^ +2 ^-^ + ^ =0 ' (A23) 

where, for example, the symmetry £ = t,rj = 1 brings us to obtain, from Eq. (|A23|i . the following restriction on the 
potential V 

i.e. a solution like this: V — exp (—20) , and therefore, <j> = hit. 
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Appendix B: Matter collineation for the scalar model with G(t) 

If the stress-energy tensor stand for a scalar model then it takes the following form: 

Tij = (p + p) UiUj + pgij , (Bl) 



where 



Then Eq. 



reads 



now reads 



C HO (G(t)T ij ) = 0, 



G p t 



P > (2 a 

t G J' 



,„ dV{4>) (2 G>\{1 V(ct>)\ 
We may follow different tactics. The first one consists in studying the whole equation 



dV(</>) (2 G'\(l V{cf>)\ 



The second one will consist in splitting the ODE in the following form (as in the standard case) 



'1 G_ 

t + 2G / 
dV(4>) _ (2 CP\ V{4>) 



in such a way that solving (|B3[) then we will be able to integrate (|B4 
Eq. (|B3[) has the following solution 



ty/GQ) 

In the same way Eq. (|B3I) admits the following symmetries: 



1 G' 



= 0, 



2 ( - t + ) ^ + m ~ 2 ^ = o. 



where the symmetry 



(B2) 

(B3) 
(B4) 



<f> = Cx I — %=+C 2 . (B5) 



t ' 2G; 6+ r^ + 2G"2G^J e + ^" 26i -°' 

1 G' \ 

- + - m + r, tt = 0, 



£ = i, 77 = -60 = 0o« , (B6) 



brings us to obtain the following constrain on function G(t) : 

G^_ 

~~G ' i 

whose solution is 
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(-112 (-11 

G"=^-±, (B7) 



G = G Q t k , keR. (B8) 

From 

Gp w r 2 Gc))' 2 »t~ 2 <=> k = 2b. (B9) 



Now, Eq. (TB4f yields 



-26-2 



whose integration gives 

lnV = -2(l + 6)ln* V = t 

such that Eq. (|B2j) is verified and therefore 

V = 4> a = (t- ab ) = t- 2b - 2 ^> a = hb+l). 

The main quantities behave as follows 

4> = 4> r b , G = Got 2b , v = V t-^ b+1 \ H = ht~ 1 , heR. (BIO) 

In the same way we also may study the following equation 

dV{4>) (2 G'\(l, V(</>)\ 

through the Lie group method. Eq. (|B11|) admits the following symmetries 



2 G' 



2 ( - + -q ) & + 2 V ^ - 4^0 = 0, 



(2 G'\ (2 G" G /2 \ 
W*U + ( ~ t + -q ) 6 + ( -J2 + -q - G2 ) z + ^ - 2& = 0, 



2 G'\ ,, . „„ /2 G' 



2 G'\ /2 G'\ , (2 G'\ ( 2 G" G 



-H 1 +a v ^ + 3 {l + G V + 2b {t + G V - t (-¥ + ^-G^^ v = ' (B15) 



(Bll) 



i4)^ +3 (M)^- 2 (?4)^-B + ? : -?) l ' ? = ' 

1 G'\ 

J + 2gJ^ = °- 
As above, the symmetry 

£ = i, 77 = -^ = o ^ b , (B12) 

brings us to obtain the following constrain on function G(t) : 

'2 G'\ / 2 G" G' 2 * 



/ G! • V *» ' G G 2 .<* = ' (B13) 
4V> + 2bV 4> - 2bV u 4> = 0, (B14) 

ill /-il2 
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Eq. (|B15|) may be rewritten as 



w.['?4W + (<» + ,,(i4)- t (-* + £-£'||r = u. 



while from Eq. (|B13I) we get 



G t 

From Eq. (|B14[) we get 



G" = ^--— , =^ G = G t k , keR. (B16) 



y^= 2 -±±Y± v = V ^ b+ V=Vot- 2 ^. (B17) 



b 

Notice that we have obtained the same results as in the splitting case 



The main quantities behave as follows 

4> = 4> r b , G = G Q t 2b , V = Vot~ 2(b+1 \_ H = ht~ 1 , heR. (B18) 
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